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Applications of negatively charged nitrogen-vacancy center in diamond exploit the center’s unique 
optical and spin properties, which at ambient temperature, are predominately governed by electron- 
phonon interactions. Here, we investigate these interactions at ambient and elevated temperatures 
by observing the motional narrowing of the center’s excited state spin resonances. We determine that 
the center’s Jahn-Teller dynamics are much slower than currently believed and identify the vital 
role of symmetric phonon modes. Our results have pronounced implications for center’s diverse 
applications (including quantum technology) and for understanding its fundamental properties. 

PACS numbers: 61.72.jn, 63.20.kp, 76.70.hb 


The negatively charged nitrogen-vacancy (NV) center 
is a point defect in diamond [l| that has found diverse 
applications in quantum technology. The center is em¬ 
ployed as a highly sensitive nanoscale sensor of electro¬ 
magnetic fields @-0], temperature [iSl and pressure 
14| that can operate in ambient and extreme conditions. 
Recent NV metrology proposals include gyroscopy [H- 

t and the development of hybrid [l^ and multi-mode 
sensors. In quantum information science, the NV 
center is used to realize spin registers |2Cll422ll a t room 
temperature and spin-photon entanglement [23l . |24 | at 
cryogenic temperatures. A new direction in NV quantum 
information science seeks to exploit spin-phonon coupling 
to enhance NV spin registers and develop novel quantum 
devices 25l-28|. 


The applications of the NV center are based upon its 
remarkable optical and spin properties. The center’s 
room temperature applications primarily rely upon its 
bright optical fluorescence, long-lived ground state spin 
coherence and methods of optical spin polarization and 
readout. The latter enable the optical detection of the 
center’s magnetic resonances (ODMR) and are the conse¬ 
quence of spin-dependent phonon-mediated intersystem 
crossings (ISCs) [i^, [s^. The center’s cryogenic appli¬ 
cations also employ the coherence of the center’s visi¬ 
ble zero-phonon line (ZPL). The necessity of cryogenics 
arises from the temperature dependent electron-phonon 
induced dephasing and depolarization of the ZPL 36, 
Electron-phonon coupling is also responsible for the mo¬ 
tional narrowing of the center’s excited state spin reso¬ 
nances, which determines their utility as an additional 
quantum resource for sensing and information process¬ 
ing 13, sH. Thus, a through understanding of the NV 
center’s electron-phonon interactions is important to the 
continued advancement of its applications and may be 
generalized to similar defects with emerging quantum ap¬ 
plications, such as the silicon-vacancy center in diamond 
Hlai and centers in silicon carbide [s^, . Here, we 

show that there exist several issues in the current under¬ 


standing and identify possible resolutions. 

The electronic structure of the NV center is depicted 
in Fig. [H The optical transitions of the visible ZPL oc¬ 
cur between the ground ^A 2 and excited spin triplet 
levels. The temperature dependent broadening of the 


ZPL was initially described 36[ using the widely applica¬ 


ble model of quadratic electron-phonon interactions with 
Ai phonon modes [s^. However, subsequent single cen¬ 
ter cryogenic measurements revealed that the broadening 
was more consistent with the characteristic oc T® temper¬ 
ature dependence of linear electron-phonon (Jahn-Teller) 
interactions with E phonon modes [s^. These inter¬ 
actions induce population transfer between the quasi¬ 
degenerate orbital states (|V), jV)) of the ^E level (see 
Fig. mi, which dephases the optical transitions and leads 
to the depolarization of the ZPL fluorescence 37| . Apply¬ 
ing their Jahn-Teller model, Fu et al 37| identified a fac¬ 
tor of ^ 2 inconsistency between the population transfer 
rates that describe the ZPL broadening and depolariza¬ 
tion at low temperatures. By introducing a phonon cutoff 
energy, Abtew et al [s^ attempted to extend Fu et al’s 
model to describe the ZPL broadening up to room tem¬ 
perature. In doing so, they obtained a cutoff at 50 meV 
for E phonons, which is much lower than the diamond 
Debye energy lud sal68 meV [i^ and the features of the 
NV electron-phonon spectral density extracted from the 
visible phonon sideband M |. 

At temperatures ^ 30 K, the complicated six level ^E 
fine structure (see Fig. [T| is observed via high resolu¬ 
tion optical spectroscopy [2^ . Above Ri 150 K [i^, the 
population transfer between the ^E orbital states is suffi¬ 
ciently fast to dynamically average the ^E fine structure 
so that ODMR resembles the simpler three level struc¬ 
ture of the ground state [i^. The dynamically averaged 
fine structure is temperature dependent and is described 
by the spin Hamiltonian 131 


H = D||(52 - -) - D^RiT)iS^, - S^) 


( 1 ) 


where Dii = 1.42 GHz and D± = 0.775 GHz are 
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FIG. 1. The electronic and fine structures of the NV center at 
high stress. The ^E sub-levels are labelled by their product of 
orbital (|X), |y)) and spin (|0), |5'±)) states, where the spin 
states are solutions of O- The ^E fine structure splittings 
{D± = D|| ± D±) are denoted in red. The ^A 2 sub-levels are 
denoted by their spin projection {nia = 0,±1). The optical 
transitions of the spin triplet and singlet levels are depicted 
as black solid arrows and the visible ZPL is at 1.956 eV. Blue 
dashed arrows represent the population transfers within the 
^E (rates Wj, and Wf). The black dashed arrows denote the 
allowed ISCs between the spin triplet and singlet levels. 


the spin-spin parameters, R{T) = _ 

-I- 1) is the temperature reduction factor, 
h and ks are the Planck and Boltzman constants, re¬ 
spectively, and is the strain splitting. Note that 
the negligible contribution of the A_l spin-orbit term (see 
Ref. (isl) to the fine structure is ignored here. 

The dynamical averaging is also expected to motion- 
ally narrow the ^E ODMR, since the rapid population 
transfer decouples the orbit and spin degrees of freedom. 
Fuchs et al have measured the ^E spin dephasing rate 
at room temperature [3]| . They attributed the observed 
dephasing to the dynamical averaging process and, us¬ 
ing a motional narrowing model, suggested that elevated 
temperatures or strain may decrease its rate. However, 
their proposal is yet to be tested by a systematic study 
of the motional narrowing effect. Furthermore, there is 
an inconsistency between Fuchs et al’s observations and 
the current ZPL broadening model. If the population 
transfer rate (^^ 10 THz) at room temperature is inferred 
from the ZPL width then the spin dephasing rate 
predicted by the motional narrowing model 1.2 MHz) 
is almost two orders of magnitude smaller than measured 
(^ 92 MHz). The anomalously low cutoff of Abtew et al, 
the discrepancy identified by Fu et al, the conspicuous 
absence of interactions with Ai modes, and the orders 
of magnitude larger than expected ^E spin dephasing 
rate, all indicate problems in the current ZPL broaden¬ 
ing model. 



FIG. 2. Example ODMR spectra at different temperatures 
(315 K upper, 455 K lower) and RF powers (440 mW left, 55 
mW right). The narrowing and reduced splitting of the lines 
at higher temperature as well as power broadening at higher 
RF power can be seen. The lineshape fits (solid lines) are the 
sum of two Lorentzians and a linear background. 


The optical polarization and readout of the spin triplet 
levels is a result of spin-selective ISCs with the inter¬ 
mediate ^Ai and ^E spin-singlet levels (see Fig. [T]) [l|. 
Goldman et al 2^, 3^ have developed a detailed model 
of the electron-phonon mechanisms that govern the ISC 
from ^E to ^Ai. However, in order to validate Goldman 
et al’s model and extend it to room temperature, more 
detailed and quantitative knowledge of linear E phonon 
interactions is required. This knowledge can be improved 
by extending the current measurements of the population 
transfer rates at cryogenic temperatures to room temper¬ 
ature and beyond. 

In this paper, we report observations of the ^E ODMR 
of NV centers in nanodiamond over the temperature 
range 295-550 K. We show that the ODMR is well de¬ 
scribed by a motional narrowing model and extract the 
population transfer rates. We establish that the rates 
are much slower than currently believed and do not ac¬ 
count for the observed ZPL broadening at room tempera¬ 
ture. We propose that quadratic Ai phonon interactions 
contribute significantly to the ZPL already above 30 K. 
Finally, we back up and rectify the proposals of Fuchs 
et al, resolve the inconsistencies of the ZPL broadening 
model and provide valuable insight into electron-phonon 
coupling above cryogenic temperatures. 

Our continuous wave ODMR experiments were per¬ 
formed using 532 nm laser excitation and fluorescence 
collection via an epifluorescence design. The nanodi¬ 
amonds were spin coated on a silica substrate. The 
NV spin resonances were driven by radio-frequency (RF) 
magnetic field created by a gold wire deposited onto the 
substrate. The excitation laser spot overlapped with the 
wire and the optical heating of the wire was used to con¬ 
trol the temperature of a chosen nanodiamond. See Ref. 


13l for further experimental details. On average, the nan¬ 
odiamonds had a diameter of ~ 30 nm and contained 
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~ 15 NV centers. We performed ODMR measurements 
on a total of 10 nanodiamonds. The results from one nan¬ 
odiamond are presented here and are consistent with the 
rest of the sample and, as will be explained, measure¬ 
ments in bulk diamond. Previous optical spectroscopy 
has measured the strain splitting of the nanodiamond 
to be /i^_L 4.7 meV M- This large strain splitting 

permits a simple fine structure (see Fig. [T]) and ap¬ 
plication of the motional narrowing model. Note that 
the strain splitting in previous reports [sil, [s^ have been 
much smaller. 

Examples of ODMR spectra are shown in Fig. [2j Av¬ 
eraging over the unresolved ^E hyperfine structure, the 
observed ODMR splitting is [l^ 

Aodmr = Id±R(T) + l[A^ + DlR^{T)] , (2) 

where A 40 MHz is the isotropic hyperfine parameter. 
RF-power broadening is evident in Fig. [2] Similar to the 
analysis of the ^^2 ODMR in Ref. l4J, a five-level model 
of the optical and spin dynamics yields the following ex¬ 
pressions for the ODMR linewidth Fqdmr and contrast 
OoDMR 


Fqdmr — -I- [ 1 -|- 


dTTAtPRF \ 

^OIOMR^l / 


1/2 


Dqdmr — Fqdmr 


AttkPrf 


AttkPff + 7i Fqdmr 


(3) 


where Fqqj^j^ and Fqq^^ir ^^e the homogenous and inho- 
mogenous linewidths in the absence of power broadening, 
Prf is the RF-power, k is a proportionality factor such 
that kPrf is the spin Rabi frequency, and 71 is the ef¬ 
fective spin relaxation rate. The essential difference to 
Ref. 13 is a much weaker, but more complicated depen¬ 
dence of 7 i on the laser power. At low laser powers, 
7 i kkisc/ik + O.b/cisc) ~ 22 MHz, where A: ft! 20 MHz 
[451 1 is the ^E radiative decay rate in nanodiamond of the 
same type and origin as used in this work and fcisc ft 50 
MHz is the average ^E ISC rate [2^|4^. Stress inhomo¬ 
geneity and the unresolved hyperfine structure contribute 
to 

^ QDMR- 

The homogenous linewidth Fqqj^jj^ = Foo + Fmn(T) 
is the sum of the broadening due to the ^E orbital de¬ 
cay rate Too = {k + 0.5fcisc)/7r and motional narrow¬ 
ing rMN(T)- Whilst the orbital decay rate increases at 
high temperature [ 1 , this temperature dependence 
is ignored in the following because the contribution of 
Too to the observed Fqdmr changes little, from 14 MHz 
to 17 MHz between 295 K and 500 K. The spin bath 
dephasing contribution to Too is also ignored because 
it has been assessed using the ^^2 ODMR to be neg¬ 
ligibly small (1-2 MHz). In the fast exchange approxi¬ 
mation of motional narrowing id^ . where the pop¬ 
ulation transfer rates iW^, W^) are much larger than 



FIG. 3. A: ODMR linewidth as a function of temperature at 
RF-powers of 400, 200, and 50 mW (top to bottom). Inset 
shows the weak optical power dependence of the linewidth 
at 294 K and at two RF powers: 47 mW (bottom) and 380 
mW (top). B and C show the RF-power dependence of the 
linewidth and contrast at 294 K and 100 mW optical power. 
D: ODMR splitting at different temperatures (50 mW RF- 
power). Error bars are determined by the statistics of re¬ 
peated measurements. The plotted linewidth is the average 
width of the two lines. 


the jump in the spin resonances between the ^E orbital 
states {2D±), Fmn(T) ft /3(T)27rIl^/FF^. The factor 
P{T) = jg giose to 1 above 

room temperatures. Thus, as IFj, increases with temper¬ 
ature, Fmn decreases. 

In the temperature regime ksT ^ Raman scat¬ 
tering of E phonons dominate the population transfer 
rates which read [ 2 ^ 

= B.T’ Jp ( 4 ) 

and where a;_L = hU/ksT and 

Of; is the cutoff energy for E phonons. The deforma¬ 
tion potential and Debye density of states for acoustic 
phonons have been assumed such that the correspond¬ 
ing electron-phonon spectral density is ~ 

and the constant Be = ^kiri^k%. Whilst in the sim¬ 
plest case V,E = lod, the cutoff is often considered as 
a phenomenological parameter which takes into account 
the departure from Je(w) oc oj^. In the high temperature 
regime 1 /!> a:_L, hflE/ksT applicable to our work, the in¬ 
tegral above equals i (l - (l + and 

= QT^, where Q is a constant. Hence, we obtain the 
final expression Fqq^j^ = Foo + l3iT)2'KD\/(QT'^). 

Systematic measurements of the ODMR linewidth, 
contrast and splitting at different temperatures, RF 
and laser powers are presented in Fig. [3l The weak 
optical-power dependence [inset of Fig. [3l[a)] supports 
the approximation 71 ft 22 MHz. The simultaneous 
fitting of the six data sets using the five parameters 
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yields = 33 ± 3 MHz, k r; 210 ± 40 MHz^ W"\ 

Cotmr = 16 ± 2%, Q = 0.83 ± 0.06 MHz and 
h^± = 4.6 ±0.2 meV. The values of k and h^± are in rea¬ 
sonable agreement with the parameters of the RF wire 
and previous optical spectroscopy, respectively. The fit¬ 
ting yields = 55 MHz at room temperature. The 

dephasing rate measured by Fuchs et al at room tem¬ 
perature in bulk diamond corresponds to FQ^j^pj^ ^ 29 
MHz. Taking into account that the much smaller stress 
splitting of Fuchs et al’s NV center will increase Q by 
~ 2, the two values are in agreement. Hence, we conclude 
that our nanodiamond measurements are consistent with 
bulk diamond and capture intrinsic phenomena of the NV 
center. 


The previous measurements of the ZPL width [37| are 
plotted in Fig. S] together with Wi/{2'k) calculated here 
using the value of Q that we obtained by fitting our 
motional narrowing observations (rescaled to = 0 to 
match the stress splitting in Ref. 13711 . It is evident that 
the rates are orders of magnitude too small to account 
for the ZPL width alone. We propose that the additional 
width is due to quadratic interactions with Ai modes 
that purely dephase the optical transitions 36, In 
which case, the ZPL width is 1491 


p. IFi IR.4 ^ 

Fzpl = -z -1-1- 7o, 

27r TT 


(5) 


where Wa is the additional dephasing rate and 70 is the 
approximately temperature independent contribution of 
the optical decay rate. As per a similar derivation of Wj,, 


ilA 

where Ba is a constant and Ha is the cutoff energy of Ai 
phonons. We used Eqs.(4-6) and fitted the ZPL width 
measurements to obtain Be = 1.32 Hz K“^, He = 13 ±1 
meV, Ba = 24 ± 4 p,Hz K“^, Ha = 37 ± 2 meV and 70 = 
16.2 ± 0.5 MHz (in bulk diamond). We confirmed our 
parameters {Be and H^;) of the population transfer rates 
by also fitting the polarization visibility measurements of 
Ref. [33 (see Fig. S]). Our fit of the visibility curve is 
practically indistinguishable from Ref. and our value 
of Be also agrees with the value ^ 1.6 Hz obtained 
there. Unlike Ref. [13, we use the same value of 
for ZPL and visibility fits and our fit to ZPL data better 
describes the low and room temperature regions than the 
extended Jahn-Teller model presented in Ref. Most 
importantly, the ZPL broadening is fully consistent with 
our ODMR measurements at elevated temperatures. 

The phonon cutoffs that we obtained are much lower 
than expected. We attribute this to the inadequacies of 
the acoustic approximation of the phonon spectral den¬ 
sity J{uj) w rjEOJ^ [H, and consider the cutoffs as 
phenomenological. Noting that He/fes ~ 155 K, these 
inadequacies are negligible at low temperatures, which 


N 

I 

CD 


D. 

N 



10 " 


Temperature (K) 


FIG. 4. Blue points are the ZPL width measured in Ref. 
[ 3 ^ . The black solid curve depicts the fit of our model and 
the black dashed curve is the extended Jahn-Teller model 
of Ref. [ 3 ^ . The red solid curve is the contribution of 
to the ZPL width according to (I4l5l) . The red dots show 
llA/27r derived from ODMR data alone. Inset: the ZPL po¬ 
larization visibility of two NV centers (red and blue points) 
from [^. The solid curve is our fit obtained using the model 
V = (Wt - IF 4 , ± r (1 - a) / (1 + o)) / (W; + W^ + r)^here 
a = 0.40 ± 0.02 and r = 80 MHz are defined in Ref. [^, and 
Wi and W-f are determined by our fit of the ZPL width. 


explains why they were not detected in previous cryo¬ 
genic measurements [ 2 ^ 371. Interestingly, is close 
to the calculated Jahn-Teller barrier energy ~ 10 meV 
[ 3 ^. Note that the spectral density extracted, for ex¬ 
ample, from the visible phonon sideband represents con¬ 
tributions of E and Ai phonons due to linear electron- 
phonon interactions. It is difficult to distinguished the 
effects of Al and E modes on the phonon band exper¬ 
imentally and ab initio calculations therefore appear to 
be the best avenue for future advancement to resolve the 
puzzle. 
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